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ABSTRACT 

We have studied a highly variable ultraluminous X-ray source (ULX) in the Fornax 
galaxy NGC 1365, with a series of 12 Chandra and XMM-Newton observations between 
2002 and 2006. In 2006 April, the source peaked at a luminosity « 3 x lO''^ erg 
s""'^ in the 0.3-10 keV band (similar to the maximum luminosity found by ASCA in 
1995), and declined on an e-folding timescale « 3 days. The X-ray spectrum is always 
dominated by a broad power-law- like component. When the source is seen at X-ray 
luminosities ~ 10"**^ erg s~^, an additional soft thermal component (which we interpret 
as emission from the accretion disk) contributes w 1/4 of the X-ray flux; when the 
luminosity is higher, ?» 3 x 10^° erg s""'^, the thermal component is not detected and 
must contribute < 10% of the flux. At the beginning of the decline, ionized absorption 
is detected around ~ 0.5-2 keV; it is a possible signature of a massive outflow. The 
power-law is always hard, with a photon index F « 1.7 (and even flatter at times), 
as is generally the case with bright ULXs. We speculate that this source and perhaps 
most other bright ULXs are in a high/hard state: as the accretion rate increases well 
above the Eddington limit, more and more power is extracted from the inner region 
of the inflow through non-radiative channels, and is used to power a Comptonizing 
corona, jet or wind. The observed thermal component comes from the standard outer 
disk; the transition radius between outer standard disk and Comptonizing inner region 
moves further out and to lower disk temperatures as the accretion rate increases. This 
produces the observed appearance of a large, cool disk. Based on X-ray luminosity 
and spectral arguments, we suggest that this accreting black hole has a likely mass 
~ 50-150Mq (even without accounting for possible beaming). 

Key words: X-rays: binaries — X-rays: individual (NGC 1365 XI) — black hole 
physics — accretion, accretion disks 



1 INTRODUCTION: ULX VARIABILITY 

It is still not clear how ultraluminous X-ray sources (ULXs) 
relate to the class of stellar-mass black holes (BHs) first 
identified and extensively studied in the Galaxy and Large 
Magellanic Cloud; in particular, whether and how they dif- 
fer from high-mass X-ray binaries (systems where the donor 
is an OB star). They could represent the upper end of the 
BH mass distribution, perhaps up to ~ lOOM©. Alterna- 
tively, the emission could be mildly beamed in our direc- 
tion, enhancing the apparent luminosity perhaps by a factor 
of a few. Finally, the emission could be super-Eddington 
by a factor of a few. One of the reasons the ULX classi- 
fication remains controversial is that very little is known 
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about their spectral state transitions. The identification 
and interpretation of state transitions over characteristic 
timescales of a few months/years has provided a powerful 
tool to understand the accretion processes (for example the 
disk/corona/jet interplay as a function of a variable mass 
accretion rate) in Galactic BHs over the last two decades 
(Remillard & McClintock 2006). 

Because of their distance, ULXs are too faint to be reg- 
ularly observed by the current generation of X-ray all-sky 
monitors, and require targeted observations. Few ULXs have 
been observed regularly enough to determine their long-term 
X-ray behaviour; most of them may vary by a factor of a few 
(compare for example the catalogues of Swartz et al. 2004 
and Liu & Bregman 2005) but tend to be persistently in a 
bright state. From this point of view, they are more similar 
to the (few) Galactic BHs in high-mass binary systems (i.e.. 
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accreting from an OB donor star) than to the more numer- 
ous ones in low-mass systems (accreting from a low-mass 
star). Significant variability or on/off transitions have been 
monitored in some of the fainter ULXs (Lx ~ a few 10"^^ erg 
s~^), for example in the Antennae (Fabbiano et al. 2003). 
However, spectral state transition or fiux variations by more 
than one order of magnitude are rarer among the brightest 
systems (Lx > 10''° erg s~'). Among the few notable excep- 
tions (apart from NGC 1365 XI, object of this paper) are the 
two brightest ULXs in M82 (Kaaret, Simet & Lang 2006), 
the brightest ULX in the Cartwheel (Wolter, Trinchieri & 
Colpi 2006), a source in NGC 3628 (Strickland et al. 2001), 
and a peculiar super-soft source that underwent an outburst 
in M 101 (Kong & Di Stefano 2005). The last source may or 
(more likely) may not belong to the same physical class as 
most other ULXs. 

In the few ULXs where the flux and/or spectral varia- 
tions are significant enough to be considered a state tran- 
sition, the pattern does not appear to correspond to the 
"canonical" sequence of physical states displayed by stellar- 
mass BHs: quiescence; low/hard state (hard power-law); 
thermal-dominant state (disk-blackbody with kTm ~ 1 
keV); very high state state (softer power-law plus disk- 
blackbody component) (Remillard & McClintock 2006; 
Fender, Belloni & Gallo 2004). Instead, ULX spectra in 
the 0.3-10 keV band seem to show three phenomenological 
features (Stobbart, Roberts & Wilms 2006), which may or 
may not correspond to truly different physical components: 
a) a broad, power-law-like component with photon index 
r « 1.7-2.3; b) a curvature or break above ~ 5 keV; c) an 
additional soft component (fcTin ~ 0.15-0.25 keV). Strong 
flux changes generally correspond to one or two of those 
features becoming more prominent or disappearing. For ex- 
ample, at higher luminosities, the power-law component ap- 
pears to increase more than the soft thermal component; 
the high-energy steepening is also more noticeable at higher 
luminosities (e.g., Roberts 2007; Dewangan et al. 2004; Kub- 
ota et al. 2001). It is difficult to classify this behaviour in 
terms of state transitions, because there is still no consensus 
on the physical origin of those three basic features. In the 
standard scenario, the soft component is emitted by a cool 
disk, the power-law component is the result of Comptoniza- 
tion, and its curvature/break is due to an increased optical 
depth of the Comptonizing region. Alternatively, it was sug- 
gested that the curved emission component may come from 
a hot accretion disk (e.g., Kubota et al. 2001; Roberts et 
al. 2005; Stobbart et al. 2006; Tsunoda et al. 2006). 

For our study, we have targeted the brightest ULX in 
NGC 1365 (henceforth, XI, following the ROSAT identifi- 
cation; Komossa & Schulz 1998) which was already known 
to be highly variable (a factor > 10) over a few months' 
timescale, from ROSAT and ASCA studies (Komossa & 
Schulz 1998; lyomoto et al. 1997). The source is located in 
the most luminous spiral galaxy of the Fornax Cluster, a su- 
pergiant barred spiral (type SBb(s)I: Sandage & Tammann 
1981) with a total gas mass ~ 3 x lO^^M©, a kinematic 
mass « 3.6 x lO'^M©, and a star formation rate ~ IOMq 
yr"^ (Roussel et al. 2001; Lindblad 1999). The Cepheid dis- 
tance to NGC 1365 is 19 ± 1 Mpc (Ferrarese et al. 2000); 
at that distance, 1" — 92 pc. In this paper, we will show 
that this source varies dramatically in luminosity over short 
timescales, and that such changes are associated to spectral 



changes, including the signature of an ionized outfiow. We 
will then discuss whether and how those changes may be 
consistent with a standard disk plus corona scenario, and 
what constraints we can obtain on the BH mass and accre- 
tion rate. 



2 SUMMARY OF OUR OBSERVATIONS 

We observed the central region of NGC 1365 with Chandra 
on 2002 December 24, then with XMM-Newton on five occa- 
sions between 2003 January 16 and 2004 July 24, and with 
Chandra again six times between 2006 April 10 and 2006 
April 23 (see Table 1 for details). Henceforth, for simplicity, 
we label our XMM-Newton/EPIC and Chandra/ ACIS obser- 
vations as E1-E5, and A1-A6, respectively, in chronological 
order. The main target was the active nucleus: findings for 
that study are reported elsewhere (Risaliti et al. 2005a,b; 
Risaliti 2007; Risaliti et al. 2007 in preparation). 



2.1 XMM-Newton observations 

For all the XMM-Newton observations, the European Pho- 
ton Imaging Camera (EPIC) was in full-window mode, with 
the medium filter, both for the pn and the two Metal Ox- 
ide Semi-conductor (MOS) detectors. We used the Science 
Analysis System (SAS) version 5.4.3 (release date 2003-09- 
11) for the EPIC spectral extraction and analysis, but we 
also checked that the results were unchanged when using 
the more recent SAS version 7.0.0 (2006-06-28). We selected 
only unflagged events (#XMMEA_EP for the pn and #XM- 
MEA_EM for the two MOSs), with pattern < 12 for the 
MOSs, pattern 4 for the pn. We defined good-time-interval 
files, removing intervals of high background. We extracted 
spectra and lightcurves for XI using the SAS task xmmse- 
lect: the source extraction region was a circle of radius 20", 
which includes « 72% of the energy in the pn and « 75% in 
the MOSs. The background was extracted from a circular re- 
gion of radius 60" , located at the same angular distance from 
the nucleus of the galaxy. We built response and auxiliary 
response files with the SAS tasks rmfgen and arfgen respec- 
tively. For each EPIC observation, we then coadded the pn 
and MOS spectra, with the method described in Page, Davis 
& Salvi (2003), in order to increase the signal-to-noise ratio. 
Finally, we fitted the background-subtracted spectra with 
standard models in XSPEC version 11.3.1 (Arnaud 1996). 
For the X-ray timing analysis, we used the SAS version 
7.0.0. We used Iccorr to correct the background-subtracted 
lightcurves for gaps in the good-time-interval. We then used 
standard XRONOS tasks for studying the lightcurves. 

In addition. Optical Monitor (OM) images were taken 
in the U, UVWl, UVM2 and UVW2 filter£|, during most 
of the observations (Table 1). We registered the astrometry 
of each OM image to a sample of bright stars in the USNO- 
Bl.O Catalog (Monet et al. 2003), so that the positional error 
is < 0'.'6 for each image, comparable to the uncertainty in 
the Chandra astrometry. We used standard IRAF tasks to 



^ See http: / / www.xmm.ac.uk/onlines / uhb /X MM_UHB / node66.htmI| 
for a plot of the OM filter throughput curves. 
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combine images in the same filter, from different observa- 
tions, (fn fact, the exposures in the UVW2 filters proved to 
be too short to provide useful information). 

2.2 Chandra observations 

Chandra observations were performed with the back- 
illuminated CCD of the Advanced Camera for Imaging and 
Spectroscopy (ACIS-S, Weisskopf et al. 2002). The first ob- 
servation was taken with the full frame; the others in the 
standard 1/4 subarray mode. The subarray mode enables a 
faster frame time (0.8 s instead of 3.2 s for full frame obser- 
vations), thus reducing the pile up fraction. This is never a 
relevant issue for the analysis of the ULX; it was selected to 
prevent high pile up fractions in the main target, the active 
nucleus. 

All Chandra datasets were analyzed with the Chandra 
Interactive Analysis of Observations (CIAO) software, ver- 
sion 3.3. For each observation, the source spectrum was 
extracted from a circular region of radius 2" centred on 
the ULX position. The background spectrum was extracted 
from circular regions near the ULX. The response and 
ancillary response functions were created using standard 
CIAO tasks, with the latest calibration files available. The 
background-subtracted spectra were later re-binned to a 
signal-to-noise ratio > 4 and fitted with standard models 
in XSPEC version 11.3.1 (Arnaud 1996). 



3 RESULTS OF OUR X-RAY STUDY OF XI 

We found three ULXs within a projected distance ~ 15 kpc 
from the active nucleus (Figures 1, 2, 3). We unofRcially 
designate them, for simplicity, XI, X2 and X3, in order of 
average X-ray brightness over the whole series of our ob- 
servations. XI is the ULX previously found by ROSAT and 
ASCA. A more accurate source position for XI, from the 
C/iandra image, is RA = 03''33™34".61, Dec = -36°09'36'.'6. 
All three ULXs (but most remarkably, XI and X3) vary sig- 
nificantly between different observations. In this paper, we 
discuss the main properties of XI, and leave a study of the 
other ULXs and non-nuclear X-ray sources in NGC 1365 to 
a follow-up paper. 



3.1 Long- and short-term variability 

The long-term variability of XI was already known from 
ASCA and _R OSAT observations (Figure 4). As a first step, 
we fitted all the Chandra and XMM-Newton spectra with 
absorbed power-law models (Table 2), confirming that its X- 
ray luminosity fluctuated by a factor of 10 over the few years 
of our observations. Our sequence of six Chandra observa- 
tions between 2006 April 10-23 caught the rise, peak and de- 
cline of an outburst, in which the source reached unabsorbed 
X-ray luminosities of (3.0 ±0.2) x 10*° erg s'^ in the 0.3-10 
keV band. This luminosity is similar to the peak observed 
by ASCA in 1995 January, which we have re-estimated from 
the original ASCA/SV& data as (3.0 ± 0.5) x erg s'^ in 
the same ban43. It was considered remarkable at that time 



^ There are inconsistencies between the ASCA count rates, fluxes 
and luminosities listed for XI in different sections of lyomoto et 



Table 4. Best-fitting parameters to the Chandra/ ACl^ spectrum 
of XI from obscrvaton Al. The XSPEC model is tbabSQal x 
tbabs X (diskbb+po). The quoted errors are the 90% confidence 
limit and Afn.Gal = 1-4 x 10^0 cm-^ (Dickey & Lockman 1990). 
The power-law normalization Kpo is in units of 10~^ photons 



Parameter 



Value 



N-ii (10^^ cm-2) 
fcTi„ (keV) 

^dbb 

r 

fiTpo (10-5) 



4.9 
0.17 
25.9 
1.83 



+4.1 
2.9 
+0.18 
0.06 
+211.0 
25.1 
+0.60 
0.60 



^■°-1.5 



xl 

^0*3-10 (10^^ erg s-l) 
^0°3-l0 (10^^ erg s-l) 



0.66(11.3/17) 

10 fi+21-5 

18.0+0^ 



that the source must have risen and declined by more that 
an order of magnitude in less than a few weeks. In fact, the 
2006 outburst shows that the characteristic timescale for the 
rise and decline is only a few days (Tables 1, 2). The peak of 
the outburst lasted for about a week, corresponding to the 
observations A3, A4 and A5 (2006 April 12-17). 

Such long-term changes in the X-ray luminosity do not 
seem to correspond to significant variations in the slope of 
the continuum. When the spectra are fitted with an ab- 
sorbed power-law model, the photon index remains ~ 1.7- 
1.9 over all our XMM-Newton and Chandra observations, 
consistent the value F « 1.7 that was estimated from the 
ASCA data during the 2005 January outburst (lyomoto et 
al. 1997). Such slopes are similar to what is often found 
in other bright ULXs, but are fiatter than the power-law 
slope in most Galactic BHs when they are in a high or very 
high state (Remillard & McClintock 2006). In fact, a more 
detailed spectral analysis based on two-component models 
(power-law plus disk-blackbody) reveals that the true slope 
of the power-law is in some cases even flatter, with a photon 
index ~ 1 (Section 3.2). 

Only two of our observations (E4 and E5) are long 
enough to enable an analyis of short-term variability. We 
find that the count rate during E5 is consistent with being 
constant, while some fiickering (quite common in accreting 
sources) is marginally detected in E4 (Kolmogorov-Smirnoff 
probability of constancy = 1.5 x IQ-'^). The total, soft-band 
and hard-band lightcurves for the EPIC-pn observation E4 
are shown in Figure 4. We do not find any evidence of X-ray 
hardness changes during that observation, nor any signifi- 
cant features in its power density spectrum. 



al. 1997. From the reported SIS count rate and fitted spectral 
parameters, we suggest that the unabsorbed luminosity must have 
been 3.2 x 10"° erg s-^ in the 0.3-10 keV band, less than half of 
what is reported in Section 2.2.1 of that paper. As a further check, 
we have re-extracted and re-analysed the ASCA data, confirming 
the lower value for the unabsorbed luminosity. 
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Table 1. Log of our observations of NGC 1365. Henceforth, we shall use the lables A#: (ACIS) and E# (EPIC) to identify the various 
observations, for simplicity. 



Label 


Date 


MJD (start) 


Mission 


Obs ID 


Good time interval (ks) 


Optical Monitor 




Al 


2002 Dec 24 


52632.62 


Chandra ACIS-S 


3554 


13.8 






El 


2003 Jan 16 


52655.98 


XMM-Newton EPIC 


0151370101 


19.0 (MOS), 15.3 (pn) 


UVW2 (3.5 ks) 




E2 


2003 Feb 09 


52679.86 


XMM-Newton EPIC 


0151370201 


4.4 (MOS), 2.5 (pn) 






E3 


2003 Aug 13 


52864.01 


XMM-Newton EPIC 


0151370701 


8.0 (MOS), 6.4 (pn) 


UVW2 (1.9 ks) 




E4 


2004 Jan 17 


53021.13 


XMM-Newton EPIC 


0205590301 


57.1 (MOS), 49.7 (pn) 


U (3.4 ks), UVWl (3.4 
UVM2 (3.4 ks) 


ks). 


E5 


2004 Jul 24 


53210.27 


XMM-Newton EPIC 


0205590401 


64.3 (MOS), 55.6 (pn) 


U (3.8 ks), UVWl (3.8 
UVM2 (4.5 ks) 


ks). 


A2 


2006 Apr 10 


53835.29 


Chandra ACIS-S 


6871 


13.4 






A3 


2006 Apr 12 


53837.52 


Chandra ACIS-S 


6872 


14.6 






A4 


2006 Apr 15 


53839.98 


Chandra ACIS-S 


6873 


14.6 






A5 


2006 Apr 17 


53842.79 


Chandra ACIS-S 


6868 


14.6 






A6 


2006 Apr 20 


53845.42 


Chandra ACIS-S 


6869 


15.5 






A7 


2006 Apr 23 


53848.42 


Chandra ACIS-S 


6870 


14.6 







Table 2. Instrumental count rates and emitted luminosities in the 0.3-10 keV band, over our series of observations (the EPIC count 
rate is the sum of pn and the two MOS's). Luminosities are based on the best fitting power-law models, for which we also list column 
densities, photon indices and normalizations (in units of 10"^ photons keV~^ cm~^ s"'^ at 1 keV). In addition, we assumed a Galactic 
line-of-sight column density of 1.4 X lO^'' cm"-^ for all fits. More complex fits were obtained for some of the observations: see Tables 3—7. 



Label Count rate 

10-2 ACIS-S cts s-i 10-2 EPIC cts s"! 



iVn (102 



Power-law fit parameters 
Photon index Kpo{10~^) L, 



0.3-10 (10^3 erg s-i) 



Al 


1.65 ±0.11 


El 




E2 




E3 




E4 




E5 




A2 


4.69 ±0.19 


A3 


6.20 ±0.21 


A4 


6.23 ±0.21 


A5 


5.30 ±0.19 


A6 


3.65 ±0.15 


A7 


1.32 ±0.10 





9 O+0.9 
^■■"•-l.C 


2.12; 


2.61 ± 0.14 




1.76; 


1.73 ± 0.25 


< 1.8 


1.75; 


2.61 ± 0.22 




2.27; 


6.49 ± 0.12 


'-'•^-0.2 


i.6i; 


6.66 ± 0.13 


3+°-2 

2 0+°-5 

1 p:+0.5 

1 6+°'^ 

1 l+°-'' 

n 5+° '' 


1.87; 
1.94; 
1.92; 
1.69; 
1.70; 
1.79; 
1.79; 



,+0.49 

'-0.40 


o 9+1.5 
■^■^-0.9 


74+1.0 

'■^-0.6 


.+0.16 
'-0.24 


1 9+0.3 
^■^-0.2 


0+0.4 

'^■■'-0.3 


+1.15 
'-0.37 




9 9 + 0.6 

^■^-0.2 


.+0.48 
-0.67 


i-3tj;i 


9 g+O-5 


+0.07 
-0.06 


2 6+°-2 
^■"-0.2 


8 K+0.4 


.+0.10 
-0.09 


9+0.3 
■^■^-0.3 


0+0.3 
''■■^-0.3 


+0.22 
^-0.22 


7 0+1.6 
'■•^-1.4 


1 8 4+^-5 
-^°-*-0.7 


,+0.17 

'-0.17 




29.6ti;J 


,+0.09 
-0.16 


9.711:^ 


29-3ti;o 


1+0.16 
-0.09 


8-3tl;^ 


24.8ti;4 


,+0.22 
-0.19 


5.7ii:t 


15.8+J;g 


,+0.53 
-0.43 




4 Q+0-5 



3.2 Thermal and non-thermal spectral 
components 

In all the observations with a sufficient number of counts, 
the X-ray spectrum appears dominated by a broad power- 
law-like component. In particular, the wider spectral cover- 
age of XMM-Newton clearly rules out the possibility that 
the source is in a thermal dominant state, when most of 
the emission comes from a standard multicolour disk. Most 
ULXs (especially those with X-ray luminosities > 10'"' erg 
s"^) seem to share this phenomenological property, indicat- 
ing perhaps that most of the disk photons are Compton- 
upscattered. 



The two longest XMM-Newton observations (E4 and 
E5) have a sufficiently high signal-to-noise ratio, allowing a 
more detailed spectral analysis. Firstly, in both datasets, sig- 
nificant systematic residuals at low energies are found when 
the source spectra are fitted with a simple power-law model: 
an additional soft component is roquireqj. The soft excess 



^ Simple power-law fits yield xl = 158/138 and 156/117 for E4 
and E5 respectively (Table 3), as opposed to xi = 129/136 and 
118/115 for disk-blackbody plus power-law fits. To verify that 
this is a statistically significant improvement, we have simulated 
pure power-law spectra using fakeit in XSPEC, and fitted them 
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Table 3. Summary of the goodness-of-fit (x5) of simple phenomenological models for each observation, and 0.3-10 keV emitted luminosi- 
ties corresponding to each of those fits. Luminosities are in units of 10^^ erg . The results for the two-component models (power-law 
-I- diskbb) and (power-law — diskbb) arc listed only for those spectra where they represent an improvement with respect to a simple 
power-law fit. Some spectra (e.g., E2) do not have enough counts for a meaningful two-component fit. 



Observation 


Parameter 


diskbb 


Spectral model 
power-law 


power-law + diskbb 


More complex models 


Al 


xi = 


1.06 (20.1/19) 


0.80 (15.3/19) 


0.66 (11.3/17) 


See Table 4 




^0.3-10 = 


4.3 


7.4 


18.0 




El 


xl = 


1.71 (37.7/22) 


0.92 (20.1/22) 


no improv 








2.2 


3.3 


— 




E2 


xl = 


2.22 (6.7/3) 


2.17 (6.5/3) 


no improv 






^-0.3-10 = 


1.7 


2.2 


- 




E3 


xl = 


1.14 (13.6/12) 


0.64 (7.71/12) 


no improv 






Lo.3-10 = 


1.7 


2.9 






E4 


xl = 


no fit 


1.15 (158.2/138) 


0.95 (129.1/136) 


See Table 5 




^0.3-10 = 


no fit 


8.5 


10.5 




E5 


xl = 


no fit 


1.33 (155.9/117) 


1.02 (117.7/115) 


See Table 6 




^-0.3-10 = 


no fit 


8.3 


9.2 




A2 


xl = 


1.15 (39.4/34) 


0.99 (33.7/34) 


no improv 






i'0.3-10 = 


11.7 


18.4 


- 




A3 


xl = 


0.82 (39.5/48) 


0.64 (30.6/48) 


no improv 


See Table 7 




^0.3-10 = 


18.9 


29.6 






A4 


Xiy 


0.80 (39.1/49) 


0.81 (39.5/49) 


no improv 


See Table 7 




^0.3-10 = 


20.1 


29.3 






A5 


Xiy 


0.71 (31.2/44) 


0.75 (32.9/44) 


no improv 


See Table 7 




Lo.3-10 = 


17.4 


24.8 






A6 




0.50 (15.5/31) 


0.38 (11.7/31) 


no improv 






Lo.3-10 = 


10.6 


15.8 






A7 


Xf 


1.28 (14.1/11) 


0.90 (9.88/11) 


no improv 






Lo.3-10 = 


2.8 


4.8 







contributes ~ 1/4 of the emitted luminositjQ in the 0.3-10 
keV band, but dominates below ~ 1.5 keV (Figure 6 and 
Tables 5, 6). This emission is well fitted by a standard disk- 
blackbody model, with colour temperatures ~ 0.3-0.4 keV. 
Such temperatures are intermediate between those gener- 
ally found in "typical" ULXs (Miller, Fabian & Miller 2004; 
Stobbart et al. 2006) and those observed from the disks 
of Galactic stellar-mass BHs (e.g., Remillard & McClintock 
2006). A possible soft excess is marginally found also for the 



both with power-law and disk-blackbody plus power-law models; 
we have then compared the decrease in xi with what we obtain 
from the real datasets. This method is essentially equivalent to 
an F-test (Protassov et al. 2002). 

^ Throughout this paper, for simplicity, luminosities are defined 
as 4:nd? times the respective fiuxes, assuming isotropic emission 
for both the power-law and the disk-blackbody components. More 
rigorously, disk-blackbody luminosities should be calculated as 
{2TTd?/ cosi) times the fiux, where cosi is the (unknown) viewing 
angle to the disk. The difference between the two estimates is less 
than a factor of 2 for i < 75° . 



Chandra observation Al (Table 4), but there are not enough 
counts to significantly constrain its flux and temperature. 

Secondly, we find that in E4 there is a significant steep- 
ening of the continuum slope above « 5.7 keV (Figure 6). 
Similar breaks are also found in various other bright ULXs: 
they have been interpreted as the effect of Comptonization 
in a warm (fcTe < 5 keV), moderately optically-thick (r ~ 2- 
4) corona (Stobbart et al. 2006; Goad et al. 2006; Dewangan, 
Griffiths & Rao 2006). The alternative interpretation for the 
spectral steepening is that the harder emission above 1.5 
keV is in fact from a hot (feTin ~ 2 keV) accretion disk, 
so that the steeper continuum above ~ 6 keV would repre- 
sent the Wien part of the thermal disk component, some- 
what modified (broadened) by Comptonization. Physically, 
this may be consistent with a "slim disk" model, also some- 
times used to explain some ULX spectra (e.g., Ebisawa et 
al. 2003; Mizuno et al. 2007). Although such model provides 
an equally good fit for the E4 spectrum, it is not accept- 
able for E5 nor for the brightest Chandra spectra; therefore 
we do not consider it as a viable physical scenario for this 
particular source. 
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Figure 1. True-colour image from our first Chandra/ ACIS-S ob- 
servation (2002 Dec 24); red corresponds to the 0.3—1 keV band; 
green to 1-2 keV; blue to 2-7 keV. On that occasions, at least 
three sources have emitted X-ray luminosities ^ 3 X 10'^^ erg s~^, 
and are labelled as XI (subject of this paper), X2, X3. The image 
was smoothed with a Gaussian kernel of radius l". 
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Figure 2. True-colour image from XMM-Newton/E,PlC (com- 
posite of pn and MOS data from the 2003-2004 observations); 
red corresponds to the 0.3-1 keV band; green to 1-2 keV; blue 
to 2-10 keV. X3 was on average much fainter than in the first 
Chandra observation (Figure 1). The image was smoothed with 
a Gaussian kernel of radius 6". 



Table 5. Best-fitting parameters to the combined XMM- 
Newton/EPlC spectrum of XI from observaton E4. The XSPEC 
models are: model 1 = tbabSQj^i X tbabs X (diskbb-|-po); model 
2 = tbabSQal X tbabs X (diskbb-|-bknpo) The quoted errors are 



the 90% confidence limit and A^jj Qj^j : 
& Lockman 1990). 



1.4 X 10^ 



(Dickey 



Parameter 



Model 1 value 



Model 2 value 



AfH (10^1 cm-2) 


1 0+°-'^ 

^■'-'-0.5 


7+°-'' 
'-0.5 


fcTin (keV) 


,-, r,q-|-0.03 

^■^•-'-0.05 


'-'■^'-0.07 


^dbb 


1 46+016 

^•^"-0.23 


n K7-I-3.20 

^•"'-0.52 


Fi 


1 07-1-0.14 
^•■-"-0.08 


1 iq-l-0.20 
^•^'-'-0.34 


Ehi 




c 7-1-0.8 
0- '-1.2 


F2 




-0.78 


Kpo (10-5) 


1 9+°-5 


1 5-1-0.5 
1-O-0.5 







0.95(129.1/136) 


0.92(122.8/134) 


rdbb 
^0.3-10 


(10^9 erg s-i) 


2 4+1'^ 


n q-l-2.6 
^■■^-0.6 


T tot 

^0.3-10 


(10^'' erg s-l) 




Q ci-l-1.3 
^•^-0.3 



Thirdly, we find that the true value of the power-law 
photon index may be lower (harder spectrum) when the 
thermal component is also taken into account. Simple power- 
law fits to the E4 and E5 spectra give F ^ 1.6-1.7 (Table 
2). But when the same E4 spectrum is fitted with a broken 
power-law plus disk-blackbody, the photon index is F « 1.1 
below the break, steepening to F ~ 2.5 at higher energies 
(Table 5). For E5, the photon index is F ^ 0.8 (Table 6), at 
least up to energies « 10 keV; we cannot determine where 
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Figure 3. True-colour UV image of NGC 1365 from XMM- 
Newton/Oyi (observed on 2004 Jan 17). Red corresponds to the 
U filter, green to the UVWl filter, and blue to the UVM2 filter. 
The positions of the nucleus and of the three ULX detected in 
the Chandra field (Figure 1) are overplotted as red circles, with a 
radius of 1'.'5, which is a conservative estimate of the uncertainty. 
North is up. East to the left. 
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Figure 4. Top panel: historical lightcurve of NGC 1365 XI. Un- 
absorbed fluxes for the ROSAT (green filled circles) observations 
have been estimated from the instrumental count rates published 
in the literature, assuming a power-law photon index F = 1.7 and 
Galactic absorption. We estimated the error bars by varying the 
spectral slope between 1.6 and 2.1, and increasing the column 
density to 2 X lO'^^ cm~'^. For the ASCA observations (magenta 
open circles) , we have re-extracted and analysed the data. For the 
Chandra (red squares) and XMM-Newton (blue triangles) obser- 
vations, we used the best-fitting models listed in Tables 3—7. We 
converted fluxes to luminosities assuming isotropic emission and 
a distance of 19 Mpc. Bottom panel: zoomed-in view of the lumi- 
nosity variations over the timespan of our Chandra observations. 



the expected spectral break occurred on that occasion. We 
tried to find evidence of high-energy spectral breaks also in 
our series of Chandra spectra. However, the ACTS sensitivity 
drops rapidly above 6 keV, making this task impossible for 
such short observations. 



3.3 The 2006 outburst sequence 

The most important result from the sequence of Chandra 
observations in 2006 April is that we caught the rise, peak 
and decline of an outburst, in which the source reached un- 
absorbed X-ray luminosities « 3 x 10*" erg s~^, similar to 
the values observed by ASCA in the 1995 January outburst. 
It was considered remarkable at that time that the source 
must have risen and declined by more that an order of mag- 
nitude in less than a few weeks. In fact, the 2006 outburst 
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Figure 5. XMM-Newton/'EPlC-pn lightcurves during observa- 
tion E4, in the total (0.3-8 keV), soft (0.3-1.5 keV) and hard 
(1.5-8 keV) bands. Short-term variability is marginally detected 
in the total band. 



Table 6. Best-fitting parameters to the combined XMM- 
Newton/EPlC spectrum of XI from observaton E5. The XSPEC 
model is tbabSQ^i X tbabs X (diskbb-|-po). The quoted errors are 



the 90% confidence limit and Nn q^i ■ 
& Lockman 1990). 



1.4 X 10^" cm"^ (Dickey 



Parameter 



Value 



A^H (10^^ cm-2) 
fcTin (keV) 

^dbb 

r 

Kpo (10-5) 



< 0.2 
-1-0.05 



0.41 

13+001 

+0.23 
0.39 
+0.4 
0.3 



0.82 
0.7 



^o°3-io (10^' erg s-i) 



1.02(117.7/115) 

2 9+05 

9 2+1° 
''•^-0.3 



shows that the characteristic timescale for the rise and de- 
cline is only a few days (Figure 4 and Tables 1, 2). The full 
sequence of spectra, fitted with absorbed power-law models, 
is shown in Figures 7-10. 

The outburst peaked between observations A3 (2006 
April 12) and A4 (2006 April 15). The best-fitting models 
are consistent with being the same for both spectra, within 
the errors: a simple power-law component (F = 1.9 ±0.2 for 
A3, F = 1.7±0.1 for A4) with neutral absorption, and essen- 
tially identical flux. In fact, the A5 spectrum is also essen- 
tially identical to A3 and A4 below « 0.7 keV and above « 2 
keV. The difference, however, is that the A5 spectrum has a 
statistically significant "soft deficit" at energies ~ 0.7-2 keV 
(Figures 8, 9), that cannot be modelled with additional cold 
absorption or with any kind of thermal emission. This fea- 
ture suggests the presence of additional ionized absorption 
in the A5 spectrum. 

To test this hypothesis, we fitted those three outburst 
spectra simultaneously, assuming the same primary spec- 
trum (same power-law slope and normalization) and the 
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Figure 6. Top panel: combined XMM-Newton/E,PlC spectrum 
and residuals from observation E4, fitted with an absorbed bro- 
ken power-law plus disk-blackbody component. The best-fitting 
parameters are given in Table 5. Bottom panel: the same for ob- 
servation E5; the best-fitting parameters are given in Table 6. 



same amount of neutral absorption, leaving only the amount 
of ionized absorption (absori model in XSPEC) free to vary 
between them. As expected, the ionized column densities in 
A3 and A4 is consistent with 0; we then fixed both of them 
at to obtain a better constraint on the error range for the 
other parameters; we also fixed the metal abundance to the 
solar value. The ionized column density of the A5 spectrum 
is « 10^^ cm~^, with a corresponding ionization parameter 
^ — Lyi3.6cv /nr^ 100, or L*' ^ 0.5-1, if we adopt the def- 
inition of ionization parameter used for example by Murray 
et al. (1995). We also note that the neutral absorbing col- 
umn density around the peak of the outburst (A3, A4, A5) 
appears to be higher (« 1.5-2 x 10^^ cm~^) than the value 
found at other epochs, when it is typically « 5 x 10^" cm~^. 

The A6 spectrum, taken three days after A5, already 
shows the beginning of the decline, with a reduction of the 
emitted luminosity by a factor « 2 (Figure 10, top panel). 
The spectrum is well fitted with a simple power law ab- 
sorbed by neutral gas. After three more days (observation 
A7: Figure 10, bottom panel) the luminosity had declined by 
another factor of ~ 3.5. In summary, the X-ray luminosity 
has an e-folding timescale ~ 3 days in the decline. 

We could compare this timescale with typical decline 
timescales observed at the end of Galactic stellar-mass BH 
outbursts. For example, the e-folding timescale for XTE 
J1550— 564 at the end of the 1998 outburst was 11 days 
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■ ^^^^ 
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Figure 7. Top panel: Chandra/ ACIS spectrum and residuals 
from observation A2. The best- fitting absorbed power-law model 
is plotted in green; fit parameters are listed in Table 2. For com- 
parison, the best-fit model from the outburst peak (observation 
A3) is overplotted in red. Middle panel: the same for observation 
A3 (peak of the outburst). Bottom panel: the same for observa- 
tion A4 (identical to A3). 



(Wu et al. 2002); and it was ^ 2 weeks for GRO J1655-40 
at the end of the 1996-1997 outburst (Sobczak et al. 1999). 
Those decline timescales very likely correspond to the vis- 
cous timescale on which the disk is emptied after the mass 
inflow rate from the donor star is strongly reduced; they 
depend on the size of the Roche lobe and on the effective 
viscosity. As the disk is emptied, a spectral state transition 
occurs from the high/soft to the low/hard state. However, 
the situation for NGC 1365 XI is entirely different. In our 
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Figure 8. Top panel: Chandra/ ACIS spectrum and residuals 
from observation A5, fitted with a power-law absorbed by both 
neutral and ionized gas. The best-fitting model is overplotted in 
green, and its parameters are given in Table 7. Bottom panel: 
expanded view, in linear scale, of the soft X-ray spectral data and 
model (plotted in green) for A5, compared with the (identical) 
best-fit model for the A3 and A4 spectra (plotted in red), to 
emphasize the effect of the ionized absorption in A5. 



interpretation, the X-ray spectrum is power-law dominated 
during the outburst; this is consistent with the suggestion 
that ULXs are in a state similar to the very high state of 
stellar-mass BHs (Roberts et al. 2006). Hence, we argue that 
the flare observed in the Chandra observations is probably 
comparable to the hard flares seen in the very high state (of- 
ten at the very beginning of an outburst) of stellar-mass BHs 
(for example, XTE J1550-564: Sobczak et al. 2000; GRO 
J1655-40: Sobczak et al. 1999; XTE J1859-H226: Brocksopp 
et al. 2002; 4U 1630-47: Trudolyubov, Borozdin & Pried- 
horsky 2001). Typical amplitudes of those X-ray flares are 
a factor of ~ 2-5; typical decline timescales are ~ 2-4 days. 
Physically, such flares may correspond to the formation and 
ejection of a Comptonizing region in the inner part of the 
accretion flow. This is consistent with the decline being as- 
sociated with an outflow. 

Finally, we checked whether there is any evidence of a 
soft excess in the 2006 outburst spectra, that can be at- 
tributed to disk emission. None of the individual spectra are 
improved by adding a disk-blackbody component. We find 
the tightest constraint by modelling the coadded spectrum 
of the observations A3 and A4. Any disk-blackbody com- 
ponent with inner-disk colour temperature feTin = 0.20 keV 
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Figure 9. Top panel: spectral residuals from observation A5, fit- 
ted with a power-law absorbed by neutral gas, with same fit pa- 
rameters as A3 and A4. Middle panel: spectral residuals when an 
additional ionized absorption component is added to the model. 
Bottom panel: model-independent ratio of the observed fluxes 
from A5 and A3, showing again the deficit of soft X-ray photons 
in A5. Dotted line: predicted ratio for an ionized absorber with 
and ^ = 100. Dash-dotted line: predicted ratio 



for Af, 



H 



10^2 cm-^ and 5 = 50. 



© 2006 RAS, MNRAS 000,[T]-?? 



10 R. Soria et al. 







' — 


H \ 1 \ — III 1 1 \ 1 \ — III : 

^ . 1„ J n n ^ 







1 2 

channel energy (keV) 



E 7 

o o 



1 


1 









0.5 1 2 5 10 

channel energy (keV) 

Figure 10. Top panel: Chandra/ ACIS spectrum and residuals 
from observation A6, fitted with an absorbed power-law. Bottom 
panel: the same for observation A7. For both spectra, the best- 
fitting parameters are listed in Table 2. The model corresponding 
to the peak of the outburst (A3) is plotted in red, for comparison. 



must have an emitted luminosity < 3.0 x 10^^ erg s~^ (90% 
confidence level) in the 0.3-10 keV band. If the inner-disk 
colour temperature fcTin = 0.40 keV, the upper limit to the 
emitted luminosity is 2.0 x 10^^ erg s~^ in the same band. 
Those apper limits are comparable to the inferred disk fluxes 
from E4 and E5. We cannot tell whether the disk has sig- 
nificantly decreased in flux during the outburst, but clearly 
it did not become more luminous. We argue that the disk 
emission represented ~ 1/4 of the X-ray emission when the 
total 0.3-10 keV luminosity was ~ 1 X 10"° erg s"^ and it 
represented < 10% of the X-ray emission when the source 
had brightened to « 3 x 10*° erg s~^. Thus, the brightening 
is entirely characterized by the enhancement of the power- 
law-like component. 



3.4 Evidence of outflows 

We further examined the behaviour at the peak of the 2006 
outburst, when spectral fitting suggested the onset of an out- 
flow. We directly compared the observed spectra A3 and A5, 
after rebinning them to increase the signal-to-noise. The ra- 
tio of the two spectra confirms in a model-independent way 
that a portion of the soft X-ray spectrum is absorbed in A5, 
consistent with an ionized column density « 10^^ cm~^, 
and an ionization parameter ^ ~ 100 (Figure 9, bottom 
panel). 



Those two characteristic values, similar to those often 
found in Seyfert galaxies (e.g., Crenshaw, Kraemer & George 
2003), are clearly inconsistent with a simple uniform shell 
of absorbing material around the X-ray source. The ionizing 
luminosity (above 13.6 eV) emitted by the ULX is « (5 ± 
1) X 10"° erg s~^. This would imply a characteristic size 
~ 5 X 10^° cm for the absorbing gas, much larger than the 
characteristic size of the binary system. Considering that the 
increase in the absorption column density has occurred over 
only few days, and that typical maximum radial velocities 
of radiatively-driven outfiows are ~ 10° cm s^^, we cannot 
expect the outfiow to have reached distances larger than a 
few xlO^'* cm. 

For a more accurate study, we tried using a radiatively- 
driven disk-wind model (Murray & Chiang 1995; Murray et 
al. 1995). We assume that L ~ I/Edd, in which case we expect 
a mass loss rate in the wind Mw ~ (lO^'^/uin) (M/IOOMq) g 
s~^. If the wind is launched from the inner disk, the initial 
(thermal) velocity Vh-i ~ 0.5-1 x 10^ cm s~^ and the mass 
loss rate ~ 10^° g s~^. The column density, integrated along 
a streamline, is 



TV 



(1) 



(Murray & Chiang 1995). For typical Voo ~ 10® cm s~^ and 
a column density ~ lO^'^ cm~^ as inferred from our spectral 
modelling, the wind would have to be launched from radii 
rin ~ 10^^ cm ~ lO'' gravitational radii for a lOOAf© BH. 
This is two orders of magnitude larger than expected, by 
comparison with disk winds in AGN. 

In summary, the problem we face is that if we were look- 
ing into a disk wind launched from rin < 10^° cm we would 
see a column density > 10^* cm~^ and, for L ~ 5 x 10'*° 
erg s~*, an ionization parameter a few orders of magnitude 
higher than inferred from the soft X-ray absorption, and 
also too high for the wind to be radiatively accelerated. In- 
creasing the mass outfiow rate would decrease the ionization 
state of the gas but increase the column density. An analo- 
gous problem was discussed by Schurch & Done (2006). A 
possible way to keep the ionization parameter in the wind 
sufficiently low is to assume (following Murray et al. 1995, 
Proga, Stone & Kallman 2000, Proga & Kallman 2004) that 
the outfiowing gas does not see most of the soft X-ray pho- 
tons: it is shielded by "hitchhiking gas" at the inner face of 
the wind, or by the disk atmosphere at the base of the wind. 
In particular, this is expected to occur when the mass out- 
flow rate is > Eddington accretion rate, because the outflow 
is Compton thick (Pounds et al. 2003) and the gas in the 
outer wind does not directly see the soft X-ray irradiation. 
However, we do see most of the X-ray flux from XI: the 
ionized absorber is only removing ~ 10% of the 0.3-2 keV 
flux from our line of sight (Figure 11). Qualitatively, we can 
explain this by assuming a low fllling factor for the absorber 
(< 10%). Alternatively, we can speculate that the observed 
soft X-ray flux has a direct (unabsorbed) component and 
an additional (~ 10%) component scattered into our line of 
sight by the outflowing gas; it is only the latter that bears 
the imprint of ionized absorption. 
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Table 7. Best-fitting parameters to the Chandra/ ACl^ spectra of 
XI from observation A3, A4 and A5, at the peak of the 2006 April 
outburst. The XSPEC model is tbabSQai X tbabs X absori X 
po. The quoted errors are the 90% confidence limit and AfH,Gal = 
1.4 X 10^0 cm-2 (Dickey & Lockman 1990). 
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4 THE MASS OF THE ACCRETING BH 

One of the unsolved mysteries regarding the nature of ULXs 
is the mass of the accreting BH. For NGC 1365 XI, there 
is realistically no chance to obtain a dynamical mass: the 
optical counterpart is not detected in the OM or in any 
other archival image we could find, and at 19 Mpc (distance 
modulus ~ 31.4 mag), a companion star would in any case 
be too faint for optical spectroscopy. Therefore, we need to 
rely on X-ray data for the physical identification of the X-ray 
source. Two arguments can be used, based on the Eddington 
limit or on the X-ray spectral features. 

The 2006 outburst peaked at an X-ray flux correspond- 
ing to an isotropic luminosity « 3 x 10'"' erg s-^ in the 
0.3-10 keV band, identical to the isotropic luminosity mea- 
sured by ASCA in 2005 January. The outburst sequence of 
our Chandra spectra suggests that the beginning of the de- 
cline is associated with an increase in ionized absorption 



(column density ~ 10^^ cm-^, ionization parameter ~ 100), 
usually indicative of an outflow. We speculate that such out- 
flow effectively blew away the X-ray emitting part of the 
inflow, causing the flux to drop by at least one order of 
magnitude. This behaviour is expected when the X-ray lu- 
minosity reaches or exceed the Eddington limit. This would 
suggest that the BH mass is < 2OOM0. Coincidentally, the 
break in the ULX luminosity function is also found at an 
X-ray luminosity « 3 x 10*'' erg s-^ (Swartz et al. 2004; 
Gilfanov, Grimm & Sunyaev 2004), perhaps suggesting that 
it corresponds to the Eddington limit of the most massive 
non-nuclear BHs in the local Universe. 



It is also possible that the emission is not isotropic, 
owing to mild geometrical beaming, so that the true lumi- 
nosity is reduced by a factor of a few (Fabrika & Mesch- 
eryakov 2001; King et al. 2001). Furthermore, even for a 
standard Shakura-Sunyaev disk, the bolometric luminos- 
ity exceeds the classical Eddington luminosity by a factor 
« 1 -h (3/5) In(M/MEdd) (Shakura & Sunyaev 1973; Begel- 
man. King & Pringle 2006; Poutanen et al. 2007). In fact, a 
steady solution exists for accretion rates M/MBdd, in which 
the excess accreting matter is blown away in an outflow, 
which is optically thin along the disc axis but optically thick 
along the other directions. This provides a natural expla- 
nation for the scattering and coUimation of the inner-disk 
emission along the disk axis (King & Pounds 2003; Begel- 
man et al. 2006; Poutanen et al. 2007). For a given observed 
flux, mild beaming and super-Eddington emission (Begel- 
man 2006) reduce the required BH mass, making it consis- 
tent with the stellar-mass range. 

The X-ray spectral analysis may provide an indepen- 
dent constraint to the BH mass and perhaps a better un- 
derstanding of the source properties. If the source is in a 
thermal-dominant state (that is, dominated by a standard, 
optically-thick, geometrically-thin accretion disk), a tight 
correlation is predicted (and observed from Galactic BHs) 
between fltted inner-disk radius, bolometric disk luminosity 
and peak color temperature of the disk: Ldisk ~ 47ri?f„T5'^. 
If we assume that i?in ~ i^isco = aGM/c? (innermost stable 
circular orbit: a = 6 in the Schwarzschild geometry, a = 1 
in the maximally-rotating Kerr geometry) , the BH mass can 
be inferred from the previous relation, as M ~ ^disk^in 
(Makishima et al. 2000; Ross, Fabian & Ballantyne 2002). 
In the case of NGC 1365 XI, a candidate disk component is 
signiflcantly detected at high signal-to-noise in the two long 
XMM-Newton observations E4 and E5 (Tables 5, 6). We ex- 
trapolate bolometric disk luminosities ~ 3.3 x 10^^ erg s-^ 
and ^ 3.5 x 10^^ erg s-^^, for E4 and E5, and fitted inner-disk 
radii of 1555/ cos 9 km and 685/ cos 9 km respectively. If we 
identify the apparent inner-disk radius with the innermost 
stable circular orbit, the BH mass would be > I8OM0 (from 
E4) or > SOMq (fi-om E5). 

This discrepancy between the two observations already 
hints that standard-disk relations may not be a reliable indi- 
cator for BH masses in ULXs. More generally, we argue that 
such disk relations should not be directly applied to ULXs 
(see also Soria & Kuncic 2007; Soria 2007). The reason is 
that the X-ray spectra of most ULXs (including NGC 1365 
XI) appear to be dominated by a Comptonization compo- 
nent (phenomenologically modelled by a power-law or bro- 
ken power-law) , not by thermal disk emission. The presence 
of a dominant Comptonized component implies that most of 
the accretion power has been removed from the disk via non- 
radiative processes, and transferred to a jet or a corona. This 
necessarily affects the direct disk emission, in particular from 
the inner region, decreasing its luminosity and temperature. 
A similar phenomenon is observed from the Galactic BH 
XTE J1550-564 in its very high state (Kubota & Done 2004; 
Done & Kubota 2006). A possible physical mechanism for 
a non-radiative, vertical extraction of power from the disk 
via magnetic torques was proposed by Kuncic & Bicknell 
(2004). Its effect is to fiatten the disk temperature profile at 
small radii; as a result, the peak of the disk emission occurs 
at larger radii (2> i?isco) and lower temperatures, leading to 
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an over-estimate of the BH mass. Another possibihty is that 
for very high accretion rates, the inner disk is completely 
replaced or covered by a moderately optically-thick corona, 
responsible for the power- law-like emission (Roberts 2007). 
In this case, too, the fitted inner-disk radius is not the inner- 
most stable circular orbit: it is instead the transition radius 
between outer (standard) disk and modified inner region of 
the inflow, much larger than Risco- Unfortunately, none of 
the brightest ULXs has ever been seen settling down into a 
disk-dominated state (unlike for example XTE J1550— 564). 
So, we do not have an independent measurement of their in- 
nermost stable orbit when the standard disk is not affected 
by Comptonization. 

In this scenario, we can use the fitted values of 7?in and 
Tin to estimate the BH mass only if we can quantify how 
the transition radius i?in moves outwards from the inner- 
most stable orbit, as a function of increasing M/MEdd (Soria 
2007). A possible suggestion is that such transition radius ~ 
spherization radius (Shakura & Sunyaev 1973) , in which case 
7?in ~ (M/MEdd)-Risco and Ti„ ~ (M/MEdd)"^^^r., where 
T* would be the peak colour temperature of a standard disk 
extended to Ribco- A more comprehensive discussion of the 
relation between transition radius and accretion rate is left 
to further work (Soria & Kuncic 2007b, in preparation). Here 
we mention only one very simple argument: if the emitted 
luminosity in the disk component is ~ 1// times the total 
luminosity (with / ~ 3-20 for "canonical" ULXs, Stobbart 
et al. 2006), we speculate that the standard disk radiates 
most of the accretion power for R > J?in ~ /-Risco, while the 
accretion power released at -Risco ^ R ^ fRiaco contributes 
to the Comptonized component. This provides a lower limit 
for the transition radius, because the radiative efficiency of 
the Comptonizing medium is less than that of the disk, and 
part of the accretion power from the inner region may also 
be carried away as mechanical luminosity of an outflow. For 
NGC 1365 XI, we found that about 3/4 of the X-ray lumi- 
nosity was in the power-law component during the E4 and 
E5 observations, which suggests that > 3/4 of the accre- 
tion power was extracted from the disk or inner inffow via 
non-radiative processes. (We also found that this fraction 
must have increased to > 90% during the 2006 outburst). 
Taken at face value, this suggests that the fitted disk radius 
Rin > 4J?isco. Thus, we argue that the X-ray spectral prop- 
erties are consistent with BH masses ~ 5O-1OOM0, even in 
the absence of beaming. 

Can we reduce this mass estimate even further? We 
need to consider that even the most luminous disks detected 
in the thermal dominant or very high states of stellar-mass 
BHs are always at least a factor of three fainter than the 
(isotropic) disk luminosities inferred for NGC 1365 XI and 
other bright ULXs. No Galactic BH has been seen with a 
disk luminosity ~ 3 x 10^^ erg s~^ as measured in XI. This 
suggests that the BH mass in this and other similar ULXs 
may still be a few times higher than in Galactic BHs, un- 
less ULX luminosities are systemetically enhanced by mild 
beaming while those of the brightest Galactic BHs are not. 



5 CONCLUSIONS AND SPECULATIONS 

We have investigated the intriguing behaviour of a ULX in 
the spiral galaxy NGC 1365 (in the Fornax cluster), char- 



acterized by high X-ray luminosity (peaking at ~ 3 x 10*" 
erg s~^) and order-of-magnitude variability over timescales 
of weeks and months. We studied the system with XMM- 
Newton and Chandra, at different luminosities between 2002 
and 2006. 

The X-ray spectrum of NGC 1365 XI is always domi- 
nated by a broad component well modelled as a power-law. 
The power-law photon index is F 1.6-1.9, consistent with 
the 1995 ASCA observations. In fact, when an additional 
thermal component is also taken into account, the slope may 
be even higher, F « 1, at least in the 1-5 keV range, with 
a possible break at higher energies. A soft excess consistent 
with disk-blackbody emission is significantly detected in the 
XMM-Newton spectra with higher signal-to-noise ratio. It 
contributes 1/4 of the observed X-ray flux when the ULX 
has a 0.3-10 keV luminosity ~ lO*" erg s~^, but it is not 
detected, and must contribute < 10% of the flux, when the 
ULX is near the outburst peak, at luminosities three times 
higher. This — and similar findings for other bright ULXs 
— suggests that one of the fundamental physical properties 
of this class of systems is that a large fraction of accretion 
power is extracted via non-thermal processes (for example, 
magnetic torques) rather than being directly radiated by an 
accretion disk, and that such fraction increases at higher 
accretion rates. A likely consequence is that the tempera- 
ture profile of the inner accretion disk may be fiatter than 
that of a standard disk; alternatively, the inner disk may 
be entirely replaced or covered by a Comptonizing region. 
In all cases, the direct disk emission component would ap- 
pear cooler and truncated at larger radii. This may lead to 
an overestimate of the central BH mass from simple spec- 
tral fitting (intermediate-mass BH scenario). A similar phe- 
nomenon has been seen in at least one Galactic BH, XTE 
J1550— 564 (Done & Kubota 2006), and may be triggered 
when M > Mfidd- 

It is sometimes noted that the spectral state behaviour 
of ULXs does not appear to follow the same pattern of 
canonical state transitions extensively studied in stellar- 
mass BHs (Section 1). We suggest that we can now start 
to outline a general sequence of spectral states applica- 
ble to both stellar-mass BHs and ULXs. Accreting BHs of 
any mass are in a low/hard state at low accretion rates 
{M < O.OlAfEdd), when the inner disk is truncated at large 
radii because not enough mass is fiowing inwards. They are 
found in a soft (thermal-dominant) state at intermediate 
accretion rates (O.lMEdd M < lAfEdd), when the disk 
reaches the innermost stable orbit and radiates all the ac- 
cretion power. They might switch to a high/hard state at 
higher accretion rates {M > MEdd) when most of the power 
begins to be extracted from the inner disk via non-radiative 
mechanisms (Foynting flux or mechanical power), and most 
of the emerging photons are coming from the Comptoniz- 
ing medium. For A'l S> A^Edd the Comptonizing corona may 
become optically thick or replace the disk altogether. The 
open question is why ULXs seem to spend most of their ac- 
tive time in the high/hard state, rather than in the thermal 
dominant state. The reason may be related either to their 
higher BH mass or, more likely, to a consistently higher mass 
transfer rate from their massive donor stars. 

An exceptional, serendipitous result of our Chandra se- 
quence of observations in 2006 April is that we caught the 
rise, peak and decline of a strong but short-lived outburst 
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(e- folding timescale ~ 3 d). This flare may be similar to 
those observed in the power-law-dominated very high state 
of Galactic BHs. There is evidence of increased absorption 
from ionized gas (column density ~ 10^^ cm~^, ionization 
parameter ~ 100) as the flux starts to decline. This might 
be due to the onset of an outflow as the source exceeded its 
Eddington limit. We speculate that the relatively low ion- 
ization parameter of the absorber is due to partial shielding 
of the wind from direct X-ray irradiation. 

Finally, we discussed two possible arguments, based on 
the X-ray luminosity and spectral features, to constrain the 
mass of the accreting BH. We argue that the most likely 
range (consistent with both arguments) is ~ SO-ISOM© even 
without taking into account the possibility of mild geometri- 
cal beaming or super-Eddington emission (which would fur- 
ther reduce the required BH mass). 
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